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I,  Project  Summary 

All  objectives  that  were  stated  in  the  Phase  I  proposal  were  met  or  exceeded  by  the  Phase  I 
effort.  Specifically,  the  following  work  was  completed; 

(i)  Development  of  an  operational  SS-MDO-U  system. 

(ii)  Introduction  of  FRM  technology  in  the  SS-MDO-U  system. 

(ill)  Validation  and  case  -  study  of  the  SS-MDO-U  development.  A  simple  example  with 
anticipated  results  was  utilized  initially  for  the  validation.  A  case  study  for  the 
simultaneous  design  of  two  ships  which  are  expected  to  operate  and  function  together  were 
analyzed  (Figure  1).  Multiple  disciplines  were  considered  during  the  design  of  each  ship 
(system  level  optimization)  and  the  combined  performance  of  the  two  ships  demonstrated 
the  system-of-systems  functionality.  This  is  a  feasibility  demonstration  of  how  the  SS- 
MDO-U  system  will  be  applied  for  designing  new  complex  concepts  such  as  the  seabasing, 
after  the  SS-MDO-U  capability  is  fully  developed. 

(iv)  Roadmaps  for  future  work  and  the  development  of  a  commercial  transition  plan  were 
completed. 


System  Level  Optimization 

Compatibility  of  seakeeping  characteristics  of  the  two  ships 


Ship  level 
Optimizations 


MDO-U  for  Heavy 
Lift  Ship 


MDO-U  for  Transport 
Ship 


Figure  1.  Flow  Chart  of  the  Case  Study  completed  during  the  Phase  I  effort 

The  SS-MDO-U  process  developed  by  MBS  allows  associating  performance  requirements  and 
economic  metrics  with  available  performance  assessment  simulation  tools  (high  or  low  fidelity). 
This  is  an  empowering  new  capability  because  it  organizes  the  complex  decision  making  process 
and  allows  to  identify  how  requirements  on  individual  disciplines  or  systems  interact  with  each 
other  and  how  they  impact  the  overall  cost.  If  this  is  not  done  properly  during  a  complex 
designing  process  then  the  cost  can  become  prohibiting. 

In  this  report  technical  information  is  presented  about  the  completed  tasks.  The  roadmap  for 
future  work  and  the  commercial  development  plan  are  included  in  the  submitted  Transition 
STTR  Plan.  Overall  this  development  sets  the  technical  foundation  for  developing  a 
computational  tool  for  assisting  Naval  engineers  during  the  decision  making  process  when 
designing  a  system  of  multiple  ships  which  must  interact  and  operate  jointly.  The  fleet  of  ships 
comprises  the  system-of-systems,  and  each  ship  comprises  each  system.  Multiple  disciplines  are 
considered  for  each  ship  during  the  design  process. 
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II,  Development  of  the  SS-MDO-U  System 

A  capability  was  developed  along  the  lines  of  the  flow  chart  presented  in  Figure  2.  It  is  based 
on  combining  multiple  MDO-U  blocks  together  for  performing  a  system-of-systems  level 
analysis.  The  flow  chart  for  each  MDO-U  block  is  presented  in  Figure  3;  it  optimizes  each 
separate  system  by  considering  the  multiple  disciplines  associated  with  the  performance  of  each 
system.  The  MDO-U  optimization  capability  has  been  developed  by  MBS  under  current  Phase  II 
SBIR  funding  from  NASA  Ames.  Previously  the  MDO-U  system  was  utilized  for  entry  vehicle 
design  by  combining  multiple  disciplines  (trajectory,  aerodynamics,  guidance  and  control, 
aerothermal,  and  thermal  analyses),  while  considering  uncertainties  in  the  thermal  loads  and  in 
the  properties  of  the  thermal  protection  systems.  In  this  project  the  MDO-U  capability 
comprised  the  building  block  for  generating  the  SS-MDO-U  system. 


Figure  2.  Flow  Chart  of  SS-MDO-U  process 

The  TC  method  is  employed  in  this  work  at  two  levels;  for  organizing  the  multiple 
optimizations  for  each  discipline  within  each  MDO-U  block;  and  on  a  system-of-systems  level 
for  organizing  the  optimizations  of  the  multiple  systems.  The  flow  chart  for  the  MDO-U  process 
is  presented  in  Figure  3. 


r. 


Figure  3.  Representative  flow  chart  of  the  Target  Cascading  method 
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At  the  top  level  of  the  optimization  proeess  in  the  TC  method  lays  the  optimization  statement 
whieh  defines  the  overall  objective  for  the  particular  ship  system.  The  top  level  optimization, 
tries  to  make  the  system’s  performance  represented  by  the  response  to  match  the  targeted 

performance  defined  as  by  the  user.  The  analysis  model  takes  into  account  the  variables 
describing  the  performance  of  the  system  ( v,  )  in  each  one  of  the  N  disciplines 
(Z)j,...,£)^)  which  are  considered.  The  analysis  model  also  takes  into  account  the  variables 
which  define  the  design  characteristics  of  the  system.  The  system  model  constitutes  a 
mathematical  statement  which  assesses  how  the  design  variables  and  the  performance  from  each 
discipline  are  combined  in  order  to  evaluate  the  overall  system  performance.  For  example,  the 
system  model  can  be  a  cost  model  which  assesses  the  value  of  a  system  based  on  the  overall  life 
cycle  cost.  In  this  manner,  when  the  expectations  for  the  performance  from  each  discipline  are 
defined  by  the  user,  the  corresponding  impact  on  the  overall  cost  can  be  identified.  Within  each 
iteration  of  the  top  level  design  optimization  process,  a  separate  optimization  problem  is  solved 
for  each  discipline.  The  linking  design  variables  are  passed  from  the  top  level 
optimization,  to  the  optimization  problem  at  each  discipline  level.  The  performance  of  each 
discipline  desired  by  the  top  optimization  problem  is  also  defined  as  . 

At  each  discipline  level  a  design  optimization  problem  is  solved  using  the  appropriate 
analysis  model.  Within  each  iteration  of  the  discipline  level  optimization  problem,  the  design 
variables  and  the  linking  variables  are  provided  to  the  simulation  model.  The  predicted 

performance  Pj^  is  returned  from  the  simulation  model.  Once  the  discipline  level  optimization 
analysis  is  completed  and  the  optimal  solution  is  reached,  the  updated  values  for  the  linking 
variables  /jj  and  the  updated  performance  Pj^  are  returned  to  the  top  level  optimization  problem. 

Then,  the  next  iteration  in  the  top  level  optimization  problem  is  performed.  The  mathematical 
statement  for  the  top  level  optimization  problem  is: 


\\P.-T,\\  +  s^ +sl 

subject  to:  XlK  “  II  -  ,  h^{Py,X^,l^)  =  Q  (1) 


’  Pd,  ’  ^D, 


1=1 


where  presents  the  inequality  constraints  and  presents  the  equality  constraints  of  the  top 

level  design  optimization  problem.  The  mathematical  statement  for  each  discipline  level  design 
optimization  problem  is: 


min  p  _  p^ 

If  A  '^A  II 

X  D,  5 

subject  to:  ,  /zj, )  <  0  ,  ,  X^  ,  )  =  0 


(2) 


The  overall  design  optimization  objective  is  achieved  by  letting  the  top  level  optimization 
(which  is  based  on  the  performance  of  each  discipline)  to  exchange  information  with  each 
discipline  in  a  very  systematic  manner.  Different  disciplines  can  share  the  same  variables  since 
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elements  of  the  groups  of  linking  variables  can  be  the  same  for  different  groups.  Thus,  the 

inter-relationships  between  disciplines  are  accounted,  while  retaining  a  very  structured  and 
efficient  manner  for  optimizing  the  design.  Another  major  benefit  of  this  method  is  that  the  user 
can  employ  different  optimization  schemes  for  each  discipline. 

In  order  to  develop  the  SS-MDO-U  process,  a  higher  level  driver  is  placed  on  top  of  the 
current  MDO-U  process  by  repeating  the  Target  Cascading  approach  at  a  higher  level.  Figure  2 
presents  details  for  the  new  SS-MDO-U  development  completed  in  this  project  in  a  flow  chart 
form.  Each  MDO-U  block  addresses  the  performance  of  each  ship,  and  the  higher  level 
optimization  process  interconnects  the  design  process  of  each  system  by  requiring  for  the  various 
ship  systems  to  operate  within  a  fleet  of  ships.  For  example,  the  seabasing  concept  represents  a 
situation  where  the  new  developments  can  be  utilized.  Linking  variables  and  performance 


metrics  are  passed  between  the  higher  level  and  the  systems’  level  design  processes.  At  the 

higher  system  of  systems  level,  the  objective  function,  the  design  variables,  and  the  constraints 
are  considered  along  with  augmentations  in  the  objective  function  introduced  by  the  limits  of  the 
constraints  imposed  by  the  TC  method.  The  mathematical  statement  for  the  higher  level  (system 
of  systems)  design  problem  is; 

mU  SSO  +  £mdo-u  +  ^  MDO-U  r3'] 

SSDV,P^o^,l^o,  ^  ^ 

subject  to; 

N 

g ss  (SSD  V ,  PfjQ, ,  ^  0  ,  hgs  (SSD  V ,  P^jj^ ,  )  =  0  ,  ^  \^md,  ~  ||  -  ^mdo-u  ’ 

(=1 


<  £ 


I 

MDO-U 


i=\ 

where  SSO  is  the  higher  level  objective  function,  SSDV  are  the  higher  level  design  variables, 
and  are  the  inequality  and  equality  constraints,  and  are  the  performance 
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functions  and  the  linking  variables  of  each  system  level  MDO-U  process  and  superscript  " 
indicates  the  values  returned  from  each  system  level  design  to  the  higher  level  system  of  systems 
design  driver. 

Uncertainty  is  introduced  in  the  SS-MDO-U  process  through  a  single  -  loop  REDO  algorithm 
for  incorporating  uncertainty  in  the  TC  method.  The  single-loop  REDO  method  imposes  the 
optimality  conditions  of  the  reliability  loops  as  an  equivalent  deterministic  equality  constraint  in 
the  design  optimization  loop.  It  therefore,  converts  the  probabilistic  optimization  formulation 
into  an  equivalent  deterministic  optimization  formulation.  It  does  not  need  to  calculate  the  MPP 
for  each  probabilistic  constraint  in  order  to  assess  the  reliability,  thus  eliminating  the  reliability 
optimization  loop  of  the  conventional  double-loop  REDO  approach.  In  general,  this  single-loop 
REDO  approach  has  comparable  efficiency  with  a  deterministic  optimization  and  the  same 
accuracy  with  the  conventional  double-loop  REDO  method.  It  is  however,  much  more  efficient 
compared  with  the  double-loop  and  decoupled  methods. 

Eased  on  the  R-percentile  formulation,  the  single-loop  REDO  problem  is  stated  by  the 
following  equivalent  deterministic  optimization  problem; 
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subject  to;  G,(d,X.,p,)  >  0  i=l,2, . n 

where 

«<  =‘'*^<^-x,p(d,X,,p,)/||<T*VG,^^p(d,X,,p,) 

d'  <  d  <  d";  Px  -  Px  -  Px 

In  Eq.  (5)  d  e  i?*"  is  the  vector  of  deterministic  design  variables,  X  e  i?'"  is  the  vector  of  random 
design  variables,  p  e  i?''  is  the  vector  of  random  and  deterministic  design  parameters  and  n,  k,  m 
and  q  are  the  number  of  constraints,  deterministic  design  variables,  random  design  variables  and 
design  parameters,  respectively.  Also,  [}^  is  the  target  reliability  index  for  the  i*  constraint, 

Px  5  Pp  are  the  mean  values  of  X  and  p,  a,  is  the  normalized  gradient  of  the  i**'  constraint  and  ct 

is  the  standard  deviation  vector  of  random  variables  X  and  parameters  p. 

The  objective  function  is  evaluated  at  the  mean  point  d,Px,Pp  and  the  constraints  are 

calculated  at  the  d,X,p  point.  The  relationship  of  Eq.  (6)  is  used  to  evaluate  the  constraints 
consistently  with  the  values  of  the  design  variables.  The  single-loop  method  does  not  search  for 
the  MPP  of  each  constraint  at  each  iteration.  Instead,  the  MPP  of  the  active  constraints  are 
correctly  identified  at  the  optimum.  This  dramatically  improves  the  efficiency  of  the  single-loop 
method  without  compromising  the  accuracy. 


(5) 

(6) 


C*10  ®  (m) 


Eigure  4.  Representative  difference  between  deterministic  and  probabilistic  optimization 

The  probabilistc  optimization  problem  stated  by  Eq.  (5)  and  (6)  is  introduced  in  the 
optimization  statements  of  the  SS-MDO-U  process  (Eigure  2).  This  has  been  accomplished  by 
subdividing  the  design  variables  and  the  linking  variables  of  the  SS-MDO-U  process  into  the 
ones  which  are  considered  deterministic  and  the  ones  which  are  considered  as  random.  The 
constraints  of  the  SS-MDO-U  process  which  contain  random  design  variables  or  random 
parameters  are  modified  according  to  Equations  (5)  and  (6)  in  order  to  account  for  the 
uncertainties.  The  modified  constraints  are  essentially  pushing  the  optimal  solution  further  back 
from  the  active  original  constraints  of  the  optimization  process.  Eigure  4  presents  graphically  the 
difference  between  a  deterministic  and  a  probabilistic  optimum.  The  probabilistic  optimal  point 
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is  placed  within  a  safety  margin  away  from  the  active  constraints  (e.g.  Point  B  vs.  Point  A  in 
Figure  4). 

During  the  Phase  I  STTR  effort  the  framework  of  the  SS-MDO-U  proeess  was  developed. 
Although  in  the  original  proposal  the  plan  was  to  just  hardwire  the  SS-MDO-U  development  in 
the  code,  a  flexible  general  purpose  SS-MDO-U  system  was  developed  instead  that  allows  the 
user  to  define  all  the  neeessary  information  for  the  system-of-systems  analysis  in  a  structured 
data  file.  This  allows  for  eonsiderably  more  flexibility  when  performing  the  ease  studies.  The 
funetionality  of  the  developed  SS-MDO-U  system  is  presented  in  Figure  5.  The  user  defines  in  a 
data  file  information  about  the  objeetive  functions,  the  linking  variables,  and  the  constraints  for 
the  system-of-systems  level,  the  various  system  levels,  and  the  various  discipline  levels 
optimizations.  The  outline  of  the  structure  of  the  data  file  is  presented  in  Figure  6. 

Operation  of  SS-MDO-U  Process 


Data  File 


SS-MDO-U  System 

Generates  optimization 
script 

automatically 

The  following  are  embedded 
in  the  optimization  script: 
•Target  Cascading  logic 
•Single-loop  RBDO 
•FRM 


Run 
optimiz^ion 
script 


SS-MDO-U 

ResiMs 


Figure  5.  Operational  flow  ehart  of  the  developed  SS-MDO-U  proeess 


The  data  files  which  are  utilized  in  the  case  study  for  the  two  ships  system  are  attached  in 
Appendix  A.  They  contain  the  definition  of  the  data  required  to  perform  the  analyses  without 
uneertainty  and  with  uncertainty,  respeetively.  During  a  potential  Phase  II  effort  the  SS-MDO-U 
eapability  ean  be  developed  further  for  becoming  more  user  friendly  and  easy  to  use.  Additional 
technieal  functionality  can  also  be  added  to  the  SS-MDO-U  system  (all  potential  new 
developments  are  outlined  in  Section  VII). 

Systom  of  Systtms  Lovol 

Design  variables,  including  any  uncertainty 

Objective  function 

Constraints 

Definition  of  each  System 

Systam  levaf  OpAnaaBon: 

Linking  variables 
Objective  function 
Constraints 
Definition  of  disciplines 

Pfsc^iBO  kmi  Op&nBatkm: 

Linking  variables 

Variables  and  parameters  with  uncertainty 
Objective  function 
Utility  functions 
Constraints 

Constraints  with  uncertainty 


Figure  6.  Strueture  of  data  file  for  SS-MDO-U  proeess 
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Ill,  Introduction  of  FRM  Technology  in  the  SS-MDO-U  System 

Capabilities  for  developing  response  surfaces  for  scalar  variables  and  time  dependent 
metamodels  for  time  dependant  responses  of  interest  are  integrated  in  the  SS-MDO-U  system. 
The  response  surfaces  for  scalar  variables  are  developed  using  the  Kriging  method  which 
comprises  a  Gaussian  type  of  metamodels.  When  the  response  of  interest  is  a  scalar  quantity  (i.e. 
roll,  heave,  or  pitch,  etc.)  then  response  surfaces  can  be  created  and  utilized  instead  of  the  actual 
simulation  models  during  the  optimization  process.  If  the  response  of  interest  is  a  time 
dependent  quantity  (i.e.  time  history  of  stress,  velocity  time  history,  etc.)  then  PCD  is  combined 
with  the  Kriging  method  for  creating  the  FRM.  The  FRM  are  employed  during  the  optimization 
process  instead  of  the  actual  solvers  for  enabling  the  use  of  high  fidelity  models  in  the  design 
process.  The  technical  background  of  the  Kriging  method  is  presented  first.  Then,  the  theory  of 
the  FRM  is  also  discussed.  Both  capabilities  have  been  integrated  in  the  SS-MDO-U  system. 

The  Kriging  metamodels  are  equations  that  replace  an  actual  simulation  model.  The  equations 
that  define  the  metamodels  are  evaluated  from  actual  simulations  performed  at  a  limited  number 
of  sample  points.  An  optimal  symmetric  latin  hypercube  algorithm  is  used  by  MBS  for  selecting 
the  sample  points.  Once  the  sample  points  and  the  values  for  the  performance  variables  are 
known,  the  Kriging  method  is  employed  for  generating  a  metamodel  for  each  performance 
function  of  interest.  The  Kriging  method  is  based  on  treating  Z(x),  the  error  between  the  actual 
performance  variable  y(x)  and  an  approximate  value  /? ,  as  a  stochastic  process; 
y  {x)=  /3+Z(x)  (7) 

where  x  is  the  d-dimensional  vector  of  the  variables  that  defines  the  point  where  the 
performance  variable  is  evaluated,  and  “d”  is  the  number  of  variables.  Z(  x  )  is  considered  as  a 
normal  process  with  zero  mean  and  a  covariance  that  can  be  expressed  as: 
cov  (Z  (x,),  Z  (xy))  =  a^R(xi,xj)  (8) 

where  cr^  is  the  process  variance  and  R(  x, ,  x^ )  is  the  spatial  correlation  function.  The  equation 

used  for  the  spatial  correlation  function  is  a  Gaussian  spatial  correlation  function; 

d  2 

)  (9) 

k=\ 

and  it  indicates  a  process  with  infinitely  differentiable  paths  in  the  mean  square  sense.  0  i-  is  the 
correlation  parameter  that  corresponds  to  the  k'*  component  of  the  d-dimensional  vector  of  the 
random  variables  x  ,  i.e.  k=l,2,. .  .,d.  For  a  set  comprised  by  “n”  number  of  sampling  points 
X  ,  1  1 ,2, . . .  ,n 

^  r=  u  .2  (10) 

The  corresponding  performance  variable  is  considered  known  and  its  values  are  defined  as; 

yl={y(x  si%y(x  ,2)v,y(^  .„)}  (ii) 

The  vector  of  correlations  between  the  sample  points  x^  and  the  evaluation  point  x  can  be 
expressed  as: 

7^(x)=  {R(y,x  ,i),R(T,T  ,2),---,  R(^  .„)}  (12) 

The  correlation  matrix  [R]  is  also  defined  among  all  the  sample  points; 

[R]  =  [R(X  ,y)]^„  (13) 

The  spatial  correlation  function  in  Equations  (12)  and  (13)  has  been  defined  by  Equation  (9).  In 
the  Kriging  method  the  value  of  the  performance  function  evaluated  by  the  metamodel  at  the 
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evaluation  point  x  is  treated  as  a  random  variable.  The  computation  of  and  Z  ( x  )  in  Equation 
(7)  is  based  on  minimizing  the  mean  square  error  (MSE)  in  the  response; 
MSE[Hv)]=E[Hv)-y(T)]2  (14) 

subjected  to  the  unbiasedness  constraint; 

E[Hx)]=E[y(x)]  (15) 

Thus,  a  set  of  parameters  0  (k=l,2,...,d)  are  computed  as  the  solution  of  minimizing  the 

product  [(det(R)) /« •  a^],  where  a^=  ^  R^\y^-F  ^),  F  is  a  vector  of  n- 

dimension  with  all  unit  entries  and  ^  defined  as; 

^={F^  R-^  F)-^  F^  R-^  (16) 


Once  the  set  of  optimal  parameters  0  has  been  computed,  then  the  performance  function  at  any 
point  x  can  be  evaluated  as; 

y{x)=p  +  7^{x)R-\y,-Fp)  (17) 


where  7^  (x)is  defined  in  Equation  (12).  Equation  (17)  comprises  the  metamodel  developed  by 
the  Kriging  method. 

The  ERM  are  developed  by  combining  PCD  with  metamodeling.  Once  they  are  created,  the 
ERM  replace  time  consuming  simulations  when  computing  the  time  history  of  various  responses 
of  interest,  and  provide  a  very  efficient  way  for  performing  the  simulations.  In  order  to  develop 
the  ERM,  the  time  histories  for  the  responses  of  interest  from  each  sample  point  are  placed  in  a 
response  matrix  [X].  The  principal  component  analysis  is  performed  for  each  response  matrix. 
Once  the  principal  components  have  been  evaluated  for  each  sample  point,  the  Kriging  method  is 
used  for  generating  the  metamodels  for  the  principal  components.  A  response  matrix  is 
comprised  by  the  time  histories  of  the  variables  of  interest; 


m= 


Xi(h)  ., 

^i(C) 

(18) 


Each  row  in  [X]  corresponds  to  a  different  location  within  the  system  which  is  analyzed,  and 
each  column  corresponds  to  a  different  time  step.  The  PCD  allows  converting  the  time  histories 
into  “modal”  type  of  information  through  singular  value  decomposition.  The  time  histories  in 
the  response  matrix  include  all  non-linear  effects.  Through  singular  value  decomposition  the 
response  matrix  [X]  is  written  as; 

(19) 

where  [U]  is  a  matrix  which  contains  “modal”  type  of  information  in  each  column.  [W]  is  a 
diagonal  matrix  which  contains  information  about  the  energy  associated  with  each  mode  placed 
in  each  column  of  [U].  [VY  is  a  matrix  which  contains  time  domain  information  about  how  the 
“modes”  are  synthesized  in  order  to  produce  the  time  histories.  The  PCD  provides  a  compact 
representation  of  the  response  for  a  non-linear  model,  and  offers  an  expansion  basis  for 
expressing  the  response  of  the  system.  In  order  to  perform  the  principal  component 
decomposition,  equation  (19)  is  re-written  as; 


[X]  =  [0 


~D 

o' 

'7' 

0 

0 

where  [W]  = 


D 

0 


0 

0 
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(20) 
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The  principal  components  decomposition  is  based  on  equation  (20)  and  results  in  the  following 
expression  for  [X]; 

[X]  =  m[D][?j]  (21) 

The  matrices  [d)^]  and  [//J  do  not  participate  in  the  definition  of  matrix  [X]  since  they  are  both 

multiplied  with  the  zero  entries  of  matrix  [W],  and  therefore  do  not  contribute  in  the  definition  of 
the  principal  components.  The  principal  components  of  the  response  matrix  are  employed  for 
developing  FRM.  Metamodels  are  developed  for  the  principal  components  of  the  response 
matrix.  The  Kriging  method  is  used  for  generating  metamodels  for  the  entries  of  the  principal 
component  matrices  resulting  into  the  equation  for  the  FRM: 

(22) 

where  y  represents  the  set  of  parameters  which  are  considered  to  vary  during  the  design  analysis. 
The  FRM  are  utilized  for  including  high  fidelity  simulations  in  the  SS-MDO-U  process. 


IV.  Simple  Validation  Example 

In  order  to  demonstrate  that  the  SS-MDO-U  system  performs  as  anticipated  a  simple  example 
is  analyzed  first.  Two  identical  single  degree  of  freedom  systems  are  analyzed.  Three 
disciplines  are  considered  for  each  single  degree  of  freedom  system;  a  static  analysis,  a  normal 
mode  dynamic  analysis,  and  a  forced  frequency  response.  Each  single  degree  of  freedom  system 
comprises  a  system  level  optimization.  At  the  system-of-systems  level  the  two  systems  are 
considered  to  be  connected  in  series  and  their  combined  static  response  comprises  the  response 
function  of  interest.  The  two  systems  are  selected  to  be  identical  at  the  starting  point  of  the 
design  in  order  to  test  that  the  final  configurations  identified  by  the  SS-MDO-U  process  are  the 
same  for  both  systems,  due  to  mathematical  symmetry. 

Figure  7  presents  the  definition  of  the  two  system  level  optimization  including  a  summary  of 
the  analyses,  the  top  level  optimization  problem,  and  the  optimization  problems  which  are 
considered  at  each  discipline  level.  It  also  depicts  the  optimization  statement  at  the  system-of- 
systems  level.  This  simple  system  is  selected  to  be  analyzed  because  its  physics  offer  a 
challenge  to  the  optimization  solution  due  to  competing  directions  and  the  limiting  constraints, 
but  yet  it  is  feasible  to  check  the  optimization  results  due  to  the  simplicity  of  the  system. 
Uncertainty  is  considered  in  the  stiffness  property  of  the  system.  This  uncertainty  impacts  the 
solution  due  to  the  constraint  which  is  imposed  on  the  maximum  allowed  natural  frequency  on 
each  one  of  the  two  systems. 

The  top  optimization  problem  at  each  system  has  as  an  objective  to  minimize  the  mass  M  of 
the  system,  while  at  the  same  time  optimizing  the  objectives  associated  with  each  one  of  the 
discipline  level  analyses.  A  lower  limit  =1  is  imposed  on  how  small  M  can  become.  The 
first  discipline  has  as  an  objective  to  minimize  the  magnitude  of  the  static  displacement  under  a 
prescribed  static  load  .  This  requirement  implies  an  expected  increase  in  the  stiffness  of  the 
system  K . 

The  second  discipline  had  as  an  objective  to  maximize  the  natural  frequency  of  the  system  (i.e. 

co^  =  )  but  with  an  upper  bound  imposed  on  the  natural  frequency  (o^  This 

optimization  problem  implies  an  expected  increase  in  K  and  a  decrease  in  M  .  However,  due  to 
the  constraint  co^  <  ,  it  is  expected  that  a  limitation  will  be  imposed  on  how  much  K  will  be 
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increased  and  how  mueh  M  will  be  redueed.  An  uneertainty  is  introdueed  in  the  value  of  the 
stiffness  K ,  whieh  is  eonsidered  as  a  random  variable  with  preseribed  standard  deviation.  Due 
to  the  variability  in  A ,  it  is  expeeted  that  the  mass  M  will  be  allowed  to  be  redueed  less  than 
before  in  order  to  still  satisfy  the  eonstraint  a>^  <  in  the  presenee  of  uneertainty.  The 


variation  in  K  will  make  a>^  take  larger  values  than  expeeted  when  K  aequires  larger  values. 

In  order  to  eompensate  for  this  it  is  expeeted  that  M  will  not  be  allowed  to  be  redueed  as  mueh 
as  before. 

The  third  optimization  diseipline  is  assoeiated  with  the  foreed  frequeney  response  analysis  of 
the  system.  A  dynamie  load  F  is  applied  on  the  system  and  the  absolute  value  of  the 


amplitude  is  minimized  (i.e. 


F 


-q}^M  +  K 


).  This  diseipline  ean  drive  the  solution  towards  two 


possible  direetions.  Either  M  will  be  redueed  and  K  will  be  inereased,  or  M  will  inerease  and 
K  will  be  redueed.  Considering  that  the  overall  objeetive  is  to  reduee  M  ,  and  that  the 
objeetives  of  the  first  two  diseiplines  are  met  only  with  an  inerease  in  K  and  a  reduetion  in  M  , 
it  is  expeeted  that  the  solution  in  the  third  optimization  diseipline  will  be  guided  towards  the 
former  path.  The  system-of-systems  level  optimization  statement  is  minimizing  the  statie 
response  of  the  system  of  the  two  springs  eonneeted  together  (i.e.  K^)  +  {H  K  2)) ) 

In  the  applieation  of  the  SS-MDO-U  eode  for  solving  this  problem  the  following  values  are 
used  for  defining  the  optimization  problem;  K^IOON  Im  and  M  =  10  kg  eomprise  the  starting 


point  for  the  optimization;  F^^  =  2,000  N  (at  both  the  system-of-systems  level,  and  at  the  two 

statie  analysis  diseiplines),  =42  is  the  limit  of  the  eonstraint,  and  F  =  200  is  the 
dynamie  exeitation.  The  results  from  the  analysis  without  eonsidering  uneertainty  and  the  ones 
with  uneertainty  are  summarized  in  Table  1.  The  results  for  both  systems  are  identieal  as 
expeeted.  In  the  analysis  with  uneertainty  the  standard  deviation  for  the  stiffness  is  eonsidered  to 
be  G  =  lOON/m  and  the  required  reliability  level  is  set  to  98%. 


initial 

value 

final  value  after  optimization 
without  uneertainty 

final  value  after  optimization 
with  uneertainty 

System-of-Systems  Level 

40 

2.269 

2.276 

System  level  (the  results  for  both  Systems  are  identieal) 

K 

100 

1,762.64 

1,757.39 

M 

10 

1 

1.113 

X, 

20 

1.134 

1.1381 

CO 

3.1623 

41.98 

39.73 

^  dyn 

0.2817 

0.1189 

0.12 

Table  1.  Summary  for  SS-MDO-U  results  from  analyses  without  and  with  uneertainty 


As  expeeted,  when  uneertainty  is  eonsidered  in  the  ealeulations,  the  amount  of  reduetion  in  the 
masses  of  the  systems  is  less  and  the  inerease  in  the  stiffnesses  is  not  as  mueh,  leading  to  a  more 
eonservative  design  that  does  not  aehieve  the  same  level  of  optimization  as  the  solution  without 
uneertainty,  but  satisfies  the  eonstraints  in  the  presenee  of  uneertainty.  When  the  stiffness 
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acquires  its  mean  value,  the  natural  frequeney  aequires  the  value  of  co^  =  39.73  Hz.  Thus,  a 
safety  margin  is  introdueed  between  the  limiting  value  of  42Hz  and  the  value  of  .  The  safety 
margin  depends  on  the  preseribed  level  of  reliability  for  meeting  the  eonstraint.  In  the  presenee 
of  uneertainty  the  solution  without  uneertainty  provides  a  natural  frequeney  equal  to  co^  =  44.37 

(i.e.  it  violates  the  eonstraint),  while  the  solution  with  uneertainty  brings  the  value  of  a>^  =  42  (i.e 
right  at  the  bounds  of  the  probabilistie  eonstraint). 


Sysi£m^/~Sj>stems  ^sponse:  Static  analysis  of  the  entire  system 
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Discipline  3:  Frequency  Response  Analysis 
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Objective:  Minimize  the  dynamic 
response  A^ofthe  system 


Figure  7.  Outline  of  the  SS-MDO-U  process  for  the  two  single  degree  of  freedom  systems 


V.  Case  Study 

A  representative  case  study  associated  with  the  simultaneous  design  of  two  ships  is  also 
analyzed  in  order  to  demonstrate  the  feasibility  of  applying  the  SS-MDO-U  development  in 
Naval  applications  of  interest  (such  as  the  Seabasing  concept).  The  case  study  accomplishes  the 
following  objectives: 

(i)  Demonstrates  how  multiple  ships  can  be  analyzed  jointly  by  considering  each  ship  as 
a  separate  system  with  multiple  disciplines  and  how  the  operational  requirements  for 
interaction  between  the  two  ships  are  considered  at  the  very  top  system-of-systems 
level. 

(ii)  Naval  architecture  disciplines  are  integrated  in  the  SS-MDO-U  process. 

(ill)  Utilizes  response  surface  Gaussian  type  metamodels  for  replacing  time  consuming 
analyses  and  for  eliminating  user  intervention  during  the  optimization  process. 

(iv)  Demonstrates  how  high  fidelity  simulations  (i.e.  UNDEX)  can  be  replaced  by  FRM 
within  the  SS-MDO-U  environment. 

(v)  Identifies  how  uncertainty  can  influence  the  solution. 
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(vi)  Demonstrates  how  eomplex  interaetions  among  diseiplines  or  among  systems  ean  be 
taken  into  aeeount  during  the  SS-MDO-U  proeess. 

The  ease  study  is  generie  but  yet  representative  of  a  multiple  ship  design  coneept  (like  the 
seabasing),  and  demonstrates  the  feasibility  of  the  SS-MDO-U  eoneept.  Figure  1  presents  the 
overall  set  up  of  the  optimization  statement.  Detailed  information  is  provided  next. 

V.l,  System-of-Systems  Level 

A  HLS  (indieated  by  subseript  “1”)  and  a  TS  (indieated  by  subseript  “2”)  are  eonsidered  to 
operate  at  the  eonfiguration  presented  in  Figure  8.  The  heading  of  the  waves  eonsidered  during 
the  simulation  is  also  presented  in  Figure  8.  A  zero  speed  and  a  sea  state  4  are  eonsidered  in  the 
simulation. 


Figure  8.  Two  ship  system 

A  eumulative  measurable  representing  the  relative  motion  between  the  two  ships  eomprises  the 
objeetive  funetion  at  the  system-of-systems  level. 

\Heavei  |  +  \Heave2 1  +  |5i  *  Roll^  |  +  |Z.2  *  Pitch2  \  (23) 


where  B  indieates  the  beam  of  the  ship,  and  L  indieates  the  length.  All  the  design  variables  of 
the  two  system  level  optimization  statements  are  eonsidered  as  design  variables  at  the  system-of- 
systems  level  (15  design  variables).  Constraints  representative  of  the  TC  method  are 
automatieally  defined  as  part  of  the  optimization  statement  aeeording  to  Equations  (4).  The 
objeetive  funetion  of  Equation  (23)  is  also  augmented  to  inelude  the  limits  of  the  TC  eonstraints 
as  part  of  the  objeetive  funetion  aeeording  to  Equation  (3). 


V.2,  System  Level  Optimization  for  HLS 

A  eost  funetion  from  the  “Prineiples  of  Eeonomies”  notes  from  the  UM  eomprises  the  objeetive 
funetion  at  the  system  level: 


minimize 


P\ 

100 


X  0.8 


(24) 


Constraints  representative  of  the  TC  method  are  automatieally  defined  aeeording  to  Equations 
(1).  The  objeetive  funetion  of  Equation  (24)  is  augmented  to  inelude  the  limits  of  the  TC 
eonstraints  aeeording  to  Equation  (1). 
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V,2,l,  Seakeeping  Discipline  for  HLS 

The  roll  performance  of  the  HLS  comprises  the  objective  function.  Minimizing  the  roll  of  the 
HLS  assists  the  objective  function  at  the  system-of-systems  level.  Constraints  are  imposed  on 
the  heave  and  pitch  responses: 

Heavey<Q.\l,  Pitchy  <Q.\55  (25) 

All  heave,  pitch,  and  roll  responses  for  both  ships  have  been  non-dimensionalized  with  2m,  2° 
and  5“,  respectively.  The  values  for  the  1/10  percentile  of  highest  values  are  considered 
throughout  the  simulations  for  the  seakeeping  responses.  The  SPP  code  from  UM  is  utilized  for 
the  seakeeping  calculations.  The  SPP  code  is  based  on  empirical  seakeeping  models  and  does 
not  account  for  the  influence  between  the  two  ships  when  evaluating  the  seakeeping  response.  A 
flow  chart  of  the  input  and  the  output  structure  of  the  SPP  is  presented  in  Figure  9. 


Figure  9.  Input  and  output  data  structure  for  SPP  code 


Gaussian  type  metamodels  are  developed  for  each  one  of  the  three  seakeeping  responses  of 
interest  using  the  Kriging  Method  (Equations  (7)  -  (17)).  This  eliminates  the  requirement  of 
using  the  SPP  code  itself  during  the  iterative  optimization  process.  This  elimination  is  necessary 
because  the  SPP  requires  from  a  user  to  interactively  adjust  certain  parameters,  thus  making  the 
direct  utilization  of  the  SPP  impossible  during  an  automated  iterative  process.  Instead,  a  limited 
number  of  seakeeping  analyses  are  performed  off-line  at  a  limited  number  of  sample  points 
selected  by  an  optimal  symmetric  latin  hypercube  algorithm.  The  results  at  the  sample  points  are 
used  for  creating  the  response  surfaces. 

The  roll,  heave,  and  pitch  responses  at  the  sample  points  are  employed  for  constructing  a 
separate  response  surface  for  each  parameter.  The  beam,  length,  draft,  and  wave  height  comprise 
the  varying  parameters  when  developing  the  response  surfaces.  Figure  10  depicts  the  three 
response  surfaces  developed  for  the  HLS  when  the  draft  and  the  wave  height  acquire  their  mean 
values  and  only  the  beam  and  the  length  are  allowed  to  vary  (this  is  done  in  order  to  visualize  the 
response  surfaces). 
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Figure  10.  Response  surfaces  for  heave,  pitch,  and  roll  of  the  HLS  (keeping  draft  and  wave 
height  at  their  mean  values  for  plotting  purposes) 

The  beam  and  the  length  of  the  HLS  comprise  the  design  variables.  However,  the  draft  of  the 
ship  is  adjusted  to  reflect  the  changes  in  the  weight  of  the  ship  through  the  changes  imposed  by 
the  structural  strength  discipline  level  analysis,  and  by  changing  the  beam  and  the  length  of  the 
ship.  The  beam  and  the  length  are  allowed  to  vary  between  120%  and  80%  of  their  original 
values. 


V,2,2  Static  Discipline  for  HLS 

A  static  analysis  is  considered  for  the  HLS  and  the  weight  is  minimized  by  retaining  the  SM 
above  a  prescribed  value.  Figure  1 1  presents  how  the  mid-ship  section  of  a  typical  HLS  is 
idealized  for  this  analysis. 


ti 


^3 

*4 

h 

Figure  11.  Simplified  mid-ship  section  used  in  static  analysis 


The  five  thicknesses  depicted  in  Figure  1 1  comprise  the  five  design  variables  of  this  optimization 
and  they  are  allowed  to  vary  between  130%  and  70%  of  their  original  values.  The  weight  is 


15 


considered  as  a  function  of  the  five  design  variables  but  it  also  depends  on  the  beam,  length,  and 
depth  of  the  HLS.  Thus,  the  outcome  of  the  seakeeping  discipline  influenees  the  static  analysis 
discipline  sinee  the  beam  and  the  length  are  determined  from  the  seakeeping  analysis.  A 
constraint  is  imposed  on  the  SM  based  on  the  maximum  possible  dimensions  which  can  be 
aequired  by  the  ship.  The  requirement  is  for  SM  not  to  be  smaller  than  227.15. 

V,3,  System  Level  Optimization  for  the  TS 

The  system  design  for  the  TS  has  been  set  up  in  an  identieal  manner  with  the  HLS.  The  main 
differences  are  in  the  dimensions  and  the  hull  form  characteristics  of  the  two  ships  and  that  an 
UNDEX  simulation  is  performed  instead  of  a  statie  analysis.  The  latter  is  done  in  order  to 
demonstrate  the  utilization  of  FRM  instead  of  high  fidelity  simulations  in  the  SS-MDO-U 
proeess.  At  the  system  level  a  cost  function  along  with  the  limits  of  the  TC  eonstraints  comprise 
the  objective  funetion  similar  to  Equations  (1)  and  (24).  The  TC  eonstraints  are  defined 
according  to  Equations  (1). 

V,3,l  Seakeeping  Discipline  for  TS 

The  seakeeping  analysis  is  similar  to  the  one  performed  for  the  HES.  The  only  difference  is 
that  the  pitch  is  minimized  in  the  objeetive  funetion  while  constraints  are  imposed  on  heave  and 
roll; 

//eave2  < 0.343 

7?o//2<0.3  ^  ' 

The  same  non-dimensionalization  has  been  imposed  to  the  heave,  pitch,  and  roll  responses  of  the 
TS  as  in  the  HES  ease.  Response  surfaces  are  created  and  utilized  in  the  optimization  proeess 
similar  to  the  HES  seakeeping  analysis.  Figure  12  presents  the  response  surfaces  for  the  TS. 


Figure  12.  Response  surfaces  for  heave,  pitch,  and  roll  of  the  TS  (keeping  draft  and  wave  height 

at  their  mean  values  in  this  set  of  results) 
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The  beam  and  the  length  comprise  the  design  variables  but  the  draft  is  also  adjusted  (similar  to 
the  HLS  seakeeping  analysis)  since  the  weight  of  the  ship  changes  from  the  UNDEX 
consideration  and  from  the  changes  in  beam  and  length.  The  complete  response  surfaces  (with 
varying  draft  and  the  wave  height)  are  utilized  during  the  optimization  process. 

V.3.2.  UNDEX  Discipline  for  TS 

A  generic  ABAQUS  FEA  model  (Figure  13)  is  utilized  for  performing  UNDEX  simulations 
for  the  TS.  The  FEA  model  for  the  ship  is  a  simplistic  one  since  the  primary  objective  is  to 
demonstrate  how  high  fidelity  simulations  can  be  included  in  the  SS-MDO-U  process  and 
developing  a  comprehensive  FEA  ship  model  is  far  beyond  the  scope  of  this  project.  In  this 
work  the  ship  is  represented  as  a  shell  structure  with  modified  modulus  of  elasticity  in  order  to 
exhibit  similar  modal  dynamic  characteristics  with  an  actual  ship.  The  ship  structure  is  divided 
into  three  sections  and  the  density  and  modulus  of  elasticity  of  each  section  comprise  the  six 
design  variables.  Both  the  ship  structure  and  the  surrounding  fluid  are  included  in  the  FEA 
model.  Infinite  finite  elements  are  employed  for  terminating  the  fluid  elements  and  for 
introducing  the  non-reflective  boundary  conditions. 


Figure  13.  FEA  model  for  TS  utilized  in  UNDEX  simulations  using  ABAQUS 

FRM  are  developed  based  on  a  complete  set  of  simulations  performed  at  a  limited  set  of  sample 
points  using  Equations  (18)  -  (22).  A  single  ABAQUS  simulation  requires  approximately  45 
minutes.  Therefore,  performing  multiple  high  fidelity  simulations  within  an  iterative 
optimization  process  is  impractical.  Instead,  the  FRM  are  constructed  and  utilized.  The  six 
design  variables  and  the  magnitude  of  the  charge  comprise  the  varying  parameters  when 
developing  the  FRM.  The  latter  allows  introducing  the  influence  of  the  changes  in  the  beam  and 
the  length  of  the  TS  to  the  loads  applied  on  the  structure  from  the  explosion.  The  changes  in  the 
beam  and  the  length  originate  from  the  seakeeping  analysis.  FRM  are  generated  for  the  time 
histories  of  the  Von  Misses  stress  at  three  different  representative  monitoring  locations  at  the 
keel.  Figure  14  presents  a  representative  correlation  between  the  FRM  and  an  actual  ABAQUS 
simulation  for  a  combination  of  the  seven  varying  variables  that  was  not  utilized  for  generating 
the  FRM.  The  results  demonstrate  a  reliable  prediction  from  the  FRM.  It  is  particularly 
important  to  notice  that  the  ability  of  MES  to  produce  FRM  is  unique  and  tailored  to  time 
dependent  responses.  The  objective  of  this  discipline  optimization  is  to  minimize  the  weight  of 
the  ship  while  constraining  the  maximum  encountered  stresses  below  the  prescribed  limit  of 
250MPa. 
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Figure  14.  Representative  eorrelation  between  FRM  and  aetual  UNDEX  simulations  for  the  time 

history  of  the  Von  Misses  stress 

V,4,  Interactions  among  Disciplines  and  Interactions  among  Systems 

The  two  diseiplines  for  the  FILS  analysis  are  interaeting  with  eaeh  other.  The  weight  of  the 
ship  whieh  ehanges  during  the  static  optimization  influences  the  seakeeping  optimization.  At  the 
same  time  the  beam  and  the  length  of  the  ship  influence  the  results  of  the  static  analysis  since 
they  participate  in  the  definition  of  the  objective  function  of  the  static  analysis.  In  a  similar 
manner  the  two  disciplines  for  the  TS  analysis  are  interdependent.  The  weight  of  the  ship  which 
changes  in  the  UNDEX  simulation  affects  the  draft  in  the  seakeeping  analysis.  At  the  same  time 
the  beam  and  the  length  of  the  ship  which  are  determined  by  the  seakeeping  simulations  affect 
the  magnitude  of  the  load  utilized  in  the  UNDEX  simulation. 

The  two  ship  systems  are  coupled  through  the  definition  of  the  objective  function  at  the 
system-of-systems  level.  The  TC  method  through  the  additional  TC  constraints  and  through  the 
augmentations  imposed  to  the  objective  functions  organize  the  multiple  optimization  problems 
which  are  being  solved  and  co-ordinate  the  changes  introduced  in  the  design  variables  at  the 
various  levels.  In  addition,  the  structured  organization  of  the  SS-MDO-U  process  that  the  TC 
method  facilitates  allows  to  account  for  the  interaction  between  the  multiple  disciplines  in  a 
structured  manner. 

V,5,  Including  Uncertainty  in  the  SS-MDO-U  Case  Study 

In  this  case  study  two  separate  system  of  systems  analyses  are  performed.  One  without 
considering  any  uncertainty  and  a  second  one  that  considers  uncertainty  in  the  wave  height 
defined  in  the  seakeeping  analyses.  The  uncertainty  influences  the  values  of  the  heave,  pitch, 
and  roll,  and  therefore  has  the  potential  of  violating  certain  constraints  in  the  seakeeping 
disciplines.  By  making  the  constraints  probabilistic  instead  of  deterministic  using  Equations  (6) 
ensures  that  the  constraints  are  satisfied  even  in  the  presence  of  uncertainty.  The  seakeeping 
constraints  become  more  conservative  when  the  single  loop  RBDO  algorithm  is  activated  and 
uncertainty  is  considered  in  the  SS-MDO-U  analysis.  The  wave  height  comprises  one  of  the 
variables  when  the  response  surfaces  are  created  for  the  heave,  roll,  and  pitch.  In  order  to 
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demonstrate  the  effect  of  the  wave  height  to  the  seakeeping  responses  the  response  surfaces  for 
the  pitch  of  the  HLS  are  presented  in  Figure  15  for  three  different  wave  heights. 
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Figure  15.  Response  surfaces  for  the  pitch  response  of  the  FILS  for  three  different  wave  heights 
(80%  of  original  value,  original  value,  and  120%  of  the  original  value,  respectively) 

V,6,  Summary  of  Results  from  Case  Study 

Table  2  summarizes  the  results  from  the  two  analyses  without  and  with  uncertainty.  The 
values  for  the  constraints  within  the  parentheses  contain  the  values  that  the  constraints  acquire  in 
the  presence  of  uncertainty.  Subscript  “1”  indicates  the  HLS  and  subscript  “2”  indicates  the  TS. 
The  thickness  variables  for  the  HLS  are  depicted  in  Figure  1 1 .  Subscripts  “1”  -  “3”  in  the 
modulus  of  elasticity  E  and  the  density  p  of  the  TS  represent  the  three  different  sections  that  the 
TS  has  been  divided  in.  All  heave,  pitch,  and  roll  responses  for  both  ships  have  been  non- 
dimensionalized  with  2m,  2°  and  5”,  respectively.  The  values  for  the  1/10  percentile  of  highest 
values  are  considered  throughout  the  simulations  for  the  seakeeping  responses  and  for  any 
objective  functions  or  constraints  that  contain  information  from  the  seakeeping  analyses.  All 
non-dimensional  variables  are  indicated  as  such  in  Table  2. 
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Initial 

configuration 

SS-MDO-U  optimum 
without  uncertainty 

SS-MDO-U  optimum 
with  uncertainty 

50m 

0.8  ND  (w.r.t.  Ri ) 

0.801  ND  (w.r.t.  Ri) 

Ly 

180m 

0.9783  ND  (w.r.t.  Li) 

1.003  ND  (w.r.t.  Li) 

h 

0.3875m 

1.1713  ND  (w.r.t.  t^) 

1.1713  ND  (w.r.t.  tj) 

h 

0.2825m 

1.1713  ND(w.r.t.t2) 

1.1713  ND  (w.r.t.  U) 

h 

0.2805m 

0.7  ND  (w.r.t.  tj) 

0.7  ND  (w.r.t.  ^3 ) 

U 

0.28m 

0.7  ND  (w.r.t.  ^4 ) 

0.7  ND  (w.r.t.  ^4  ) 

h 

0.335m 

0.7  ND  (w.r.t.  tj) 

0.7  ND  (w.r.t.  tj) 

B2 

25.4m 

0.8  ND  (w.r.t.  R2) 

0.8589  ND  (w.r.t.  B^) 

Li 

175m 

0.8  ND  (w.r.t.  L2  ) 

0.8  ND  (w.r.t.  L2  ) 

Ey 

3.0E+14Pa 

1.0  ND  (w.r.t.  E^) 

1.0  ND  (w.r.t.  E^) 

El 

3.0E+14Pa 

1.0  ND  (w.r.t.  E^) 

1.0  ND  (w.r.t.  £3) 

El 

3.0E+14Pa 

1.0  ND  (w.r.t.  E^) 

1.0  ND  (w.r.t.  £3) 

P\ 

7,800kg/m" 

0.702  ND  (w.r.t.  ) 

0.702  ND  (w.r.t.  p^ ) 

Pi 

7,800kg/m^ 

0.7  ND  (w.r.t.  P2) 

0.7  ND  (w.r.t.  P2 ) 

Pi 

7,800kg/m" 

0.7  ND  (w.r.t.  P2 ) 

0.7  ND  (w.r.t.  P2 ) 

System-of  Systems  Objective  Function, 

ND 

1 

0.7417 

0.7677 

FiLS  system  level  Objective  Function,  ND 
Ship  Cost 

1 

0.8221 

0.8387 

HLS  Seakeeping  Discipline  Objective 
Function,  Roll,  ND 

0.382 

0.2057 

0.2113 

HLS  Seakeeping  Discipline,  Constraint 
on  Heave  <0.17,  ND 

0.132 

0.1458(0.1524) 

0.1545  (0.1588) 

HLS  Seakeeping  Discipline, 

Constraint  on  Pitch  <0.155,  ND 

0.117 

0.1399  (0.1815) 

0.1118(0.155) 

HLS  Static  Discipline,  Objective  function. 
Weight,  ND 

1 

0.826 

0.8469 

HLS  Static  Discipline,  Constraint  on 

SM<  227.15 

252 

227.8 

228.2 

TS  System  Level  Objective  Function,  ND 
Ship  Cost 

1 

0.6997 

0.7407 

TS  Seakeeping  Discipline  Objective 
Function,  Pitch,  ND 

0.317 

0.3009 

0.3147 

TS  Seakeeping  Discipline,  Constraint  on 
Heave  <  0.343  ,ND 

0.3054 

0.265  (0.3075) 

0.2638  (0.306) 

TS  Seakeeping  Discipline,  Constraint  on 
Roll  <  0.3,  ND 

0.2274 

0.184  (0.2752) 

0.1702  (0.2809) 

TS  UNDEX  Discipline,  Objective 

Function,  Weight,  ND 

1 

0.4484 

0.4810 

TS  UNDEX  Discipline,  Maximum  Von 
Misses  stress  in  Section  1  <  250 

108.68 

72.6 

76.76 

TS  UNDEX  Discipline,  Maximum  Von 
Misses  stress  in  Section  2  <  250 

207.13 

104.74 

112.69 

TS  UNDEX  Discipline,  Maximum  Von 
Misses  stress  in  Section  3  <  250 

83.37 

37.84 

40.606 

Table  2.  Summary  of  Results  from  Case  Study;  SS-MDO-U  analyses  of  two  ships;  without 

uneertainty  and  with  uneertainty 
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The  main  observations  from  the  two  SS-MDO-U  analyses  are: 

(i)  There  is  a  signifieant  improvement  (-25%)  in  the  system-of-systems  objeetive  funetion  that 
represents  the  maximum  relative  motion  between  the  two  ships. 

(ii)  There  is  a  signifieant  reduetion  (-20%  and  -30%)  in  the  cost  associated  with  both  ships  at 
each  one  of  the  two  system  levels  objective  functions. 

(iii)  Due  to  the  dependency  of  the  cost  functions  on  the  main  dimensions  of  the  ships,  the  beam 
and  the  length  of  both  ships  is  reduced  as  much  as  the  lower  bounds  of  the  design  variables 
and  as  much  the  seakeeping  constraints  allow.  In  addition,  the  requirement  for  minimizing 
the  relative  motion  between  the  two  ships  at  the  system-of-systems  level  influences  the 
reduction  of  the  beam  and  length  variables. 

(iv)  The  constraint  on  the  pitch  of  the  HLS  is  violated  in  the  presence  of  uncertainty  from  the 
deterministic  optimum  configuration.  When  uncertainty  is  included  in  the  SS-MDO-U 
process  the  optimal  solution  is  adjusted  and  the  constraint  is  satisfied  even  in  the  presence 
of  uncertainty. 

(v)  The  variables  for  the  modulus  of  elasticity  of  the  TS  do  not  change  from  their  original 
values  in  the  UNDEX  discipline.  That  makes  sense  because  the  constraints  (which  are 
influenced  by  the  E)  do  not  become  active  and  the  objective  function  does  not  depend  on  E, 
therefore  there  is  no  reason  for  inducing  any  changes.  The  densities  are  associated  with  the 
weight  of  the  ship  and  they  are  pushed  all  the  way  to  their  lower  bound  since  the  stress 
constraints  remain  inactive.  The  reason  for  the  latter  is  the  reduction  in  the  total  pressure 
load  from  the  UNDEX  due  to  the  reduced  footprint  of  the  TS  which  is  created  from  the 
reduction  in  beam  and  length. 

Overall  the  SS-MDO-U  process  operates  as  expected.  It  leads  to  improved  designs  at  the 
system-of-systems  level  and  at  each  system  level.  It  also  accounts  properly  for  the  interactions 
among  the  multiple  systems  and  among  the  multiple  disciplines  within  each  system. 

Eigure  16  presents  on  the  response  surfaces  of  the  heave,  pitch,  and  roll  of  the  HES  (under 
constant  draft  and  wave  height  for  visualization  purposes,  in  reality  these  two  variables  change 
as  well)  the  starting  point,  an  intermediate  point,  and  the  final  points  for  the  optimized 
configurations  without  and  with  uncertainty.  In  the  response  surface  for  the  pitch  it  can  be 
observed  how  the  optimal  point  under  uncertainty  has  been  adjusted  to  a  lower  pitch  response  in 
order  to  satisfy  the  constraint  in  the  presence  of  uncertainty.  The  two  other  responses  increase 
by  small  amounts  in  order  to  accommodate  this  change,  but  still  remain  within  the  feasible  limits. 
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Figure  16.  Tracing  the  starting  Point  “0”,  intermediate  point  “3”,  deterministic  optimum  “8”,  and 
probabilistic  optimum  “U”  in  the  seakeeping  response  surfaces  of  the  HLS  (the  response  surfaces 
are  presented  under  constant  draft  and  wave  height  for  visualization  purposes,  during  the 
optimization  these  two  variables  vary  as  well) 

VI,  Contributions  of  this  Project 

The  following  are  the  new  contributions  of  the  completed  Phase  1  project: 

•  Developed  a  general  purpose  SS-MDO-U  code. 

•  Introduced  response  surface  and  FRM  capabilities  in  the  SS-MDO-U  development. 

•  Demonstrated  through  a  simple  analysis  that  the  SS-MDO-U  process  operates  as  expected. 

•  Demonstrated  through  a  generic  but  relevant  to  Naval  interests  case  study  that: 

(i)  Interactions  among  multiple  ship  systems  can  be  accounted  properly. 

(ii)  Naval  Architecture  disciplines  can  be  included  in  the  SS-MDO-U. 

(iii)  High  fidelity  simulations  can  be  utilized  in  the  SS-MDO-U  process. 

(iv)  Interactions  among  disciplines  for  a  particular  ship  are  accounted  properly  within 
each  system  analysis. 

(v)  Interaction  and  compatibility  among  systems  is  coordinated  at  the  system-of-systems 
level. 

(vi)  Physical  attributes  driven  by  performance  requirements  can  be  linked  to  cost 
objectives  within  the  SS-MDO-U  process,  thus  optimizing  performance  while 
reducing  the  cost. 
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(vii)  Accounting  for  uncertainties  during  the  SS-MDO-U  simulations  introduees 

eonservatism  in  the  eonstraints  and  ensures  that  the  optimal  point  will  satisfy  the 
eonstraints  even  in  the  presenee  of  uneertainty. 

The  TC  approaeh  dietates  the  strueture  of  the  SS-MDO-U  system.  Overall,  the  TC  divides  a 
very  eomplex  and  eomplieated  optimization  ehallenge  into  a  group  of  smaller  and  systematieally 
intereonneeted  optimization  analyses.  It  allows  to  proeess  the  results  from  eaeh  system  and  eaeh 
diseipline  separately  and  thus  identifies  the  eritieal  faetors  in  the  deeision  making  proeess.  The 
TC  approaeh  provides: 

(i)  a  symmetrie  organization  to  the  diseipline  level  optimizations  within  eaeh  system 
analysis 

(ii)  a  systematie  flow  of  data  among  the  diseiplines  of  eaeh  system 

(iii)  a  system  of  organized  links  between  multiple  systems  at  the  system-of-systems  level 
analysis 

(iv)  a  eapability  to  drive  the  optimal  point  to  a  more  eonservative  solution  in  order  to 
avoid  violation  of  all  eonstraints  in  the  presence  of  uneertainty 


VII.  Potential  New  Developments  and  Market  Potential 

The  Phase  I  STTR  effort  established  the  eore  competeney  of  the  SS-MDO-U  eapability.  It 
also  demonstrated  the  feasibility  of  the  new  teehnology  through  a  ease  study.  The  eore  value  of 
the  SS-MDO-U  is  that  it  aeeounts  for  interaetions  and  eommon  physieal  attributes  between 
multiple  diseiplines  within  a  system,  it  integrates  multiple  systems  for  a  system-of-systems  level 
analysis  and  it  relates  physieal  eharaeteristies  to  eost.  It  provides  the  enabling  teehnology  for 
dividing  a  very  eomplex  and  intereonneeted  design  problem  into  a  network  of  multiple  smaller 
and  systematieally  linked  optimization  problems,  thus  making  the  solution  feasible.  At  the  same 
time  it  allows  to  identify  the  eritieal  points  and  parameters  to  the  overall  deeision  making  proeess. 
During  the  Phase  II  effort  the  following  Tasks  will  be  pursued: 

(i)  Utilize  game  theory  in  order  to  eonvert  desired  operational  and  performanee  eharaeteristies 
into  mathematieal  statements  that  ean  be  employed  in  the  mathematieal  artieulation  of  an 
optimization  statement. 

(ii)  Develop  an  algorithm  whieh  identifies  the  feasibility  of  the  starting  seleetion  with  respeet 
to  all  of  the  eonstraints.  Sinee  multiple  designs  and  multiple  diseiplines  within  eaeh  design 
are  eonsidered  during  the  system  of  systems  optimization  proeess,  it  is  important  to  ensure 
that  the  initial  starting  seleetion  of  design  variables  resides  within  the  feasible  region  for  all 
eonstraints  from  all  diseiplines  and  all  individual  ship  platforms. 

(iii)  Develop  an  algorithm  that  adjusts  the  final  optimum  in  order  to  ensure  that  all  the 
eonstraints  from  the  lower  level  diseiplines  and  systems  are  satisfied.  Sinee  optimizations 
are  performed  at  multiple  levels  of  the  SS-MDO-U  proeess,  the  adjustment  of  the  design 
variables  is  performed  at  multiple  levels.  The  development  under  this  task  will  ensure  that 
all  of  the  eonstraints  of  the  various  level  optimizations  will  be  satisfied  at  the  end  of  the  SS- 
MDO-U  proeess. 

(iv)  Identify  the  most  important  variables  whieh  have  been  eneountered  during  the  optimization 
proeess.  This  eapability  will  provide  information  to  the  user  about  the  nature  of  the 
optimization  that  is  been  solved  and  the  relative  importanee  of  the  various  variables  utilized 
in  the  proeess. 
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(v)  Create  a  user  interfaee  whieh  will  guide  the  user  in  ereating  the  data  file  for  the 
optimization  proeess  and  faeilitate  the  visualization  of  the  results.  In  order  to  inerease  the 
utilization  of  the  new  produet  a  “wizard”  based  interfaee  will  guide  the  user  in  the  required 
definition  of  the  design  variables,  objeetive  funetions,  and  eonstraints  required  for  the 
mathematieal  artieulation  of  the  SS-MDO-U  optimization  proeess. 

(vi)  Enhanee  the  library  of  optimization  algorithms  whieh  ean  be  employed  for  guiding  the 
multiple  optimization  proeesses.  This  will  provide  additional  flexibility  in  definition  of  the 
optimization  statements. 

(vii)  Develop  a  eapability  for  handling  diserete  design  variables. 

(viii)  Perform  a  ease  study  representative  of  the  Seabasing  eoneept. 

The  market  potential  is  extensively  discussed  in  the  submitted  STTR  Transition  Plan. 
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Appendix  A. 

SS-MDO-U  data  files  for  the  case  study  analyses  without 
uncertainty  and  with  uncertainty 
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ca5e_5t  udy[ 1] ,  t  xt 

$SS-MDO-U  Analysis  without  Uncertainty 
$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

$  Define  system  of  systems  level  optimization 
BEGI  N  SOS 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  system  of  systems  level  optimization  -  SOS_Main,m 
BEGI  N  MAI  N 

$ 

=^DESI  GNVARI  ABLE 

$  name  initial-value  lower-bound  upper-bound  nondimensionalization-facor 
Standard-devi  ati  on 


Bl 

50,  0 

0 

8H0,  0 

1 

2H0,  0 

50,  0 

LI 

180,  0 

0 

8H80,  0 

1 

2H80,  0 

180,  0 

t  1 

0,  3875 

0 

7H,  3  8  7  5 

1 

3H,  3  8  7  5 

0,  3875 

t  2 

0,  2825 

0 

7H,  2  8  2  5 

1 

3H,  2  8  2  5 

0,  2825 

t  3 

0,  2805 

0 

7H,  2  8  0  5 

1 

3H,  2  8  0  5 

0,  2805 

t  4 

0,  2800 

0 

7H,  2  8  0  0 

1 

3H,  2  8  0  0 

0,  2800 

t  5 

0,  335 

0 

7H,  3  3  5  0 

1 

3H,  3  3  5  0 

0,  3350 

B2 

25,  4 

0 

8H5,  400 

1 

2H5,  400 

25, 400 

L2 

175,0 

0 

8H75,  00 

1 

2H75,  00 

175,00 

El 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

E2 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

E3 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

r  hoi 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

r  ho2 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

rho3 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

$ 

=^OPTi  Ml  ZATi  ON  OPT!  ON 
$  opt  i  on  value 

Di  s  pi  ay  iter 

Diagnostics  off 
LargeScale  off 
MaxFunEvals  300 
DiffMinChange  0,005 
DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  systemof  systems  level  optimization  - 
SOS  Obj  ect i  ve,  m 
BEGl  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 

WH 1  2,2 

Wei  ght  0  (  50,  0/  2M  0,  3  8  7  5  +0,  2  8  2  5  )  +1  2  M  0,  2  8  0  5  +0,  2  8  0  0  +0,  3  3  5  0  )  )  M80,  0 
Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 
DispH  6  1  2  3  9,9 

Displ  DispHMWeightl-WeightO)/WeightO 

TDl  Di  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

WH2  2,2 

A  1  3  6  9  8,  8  2  H2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  1  7  H2H  2  /  2  5,  4  /  1  7  5,  0 

A3  3  8  4  8,  2  0  H2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  OH  6  9  8,  8  2  +7  8  0  0,  0  H  4  4  9,  1  7  +7  8  0  0,  OH  8  4  8,  2  0 

Massl  r  h  0 1  HI +r  h  0  2  H2 +r  h  0  3  H3 
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Di  spT  2  2  9  2  9,  1 

Di5p2  Di  5  pTM  Ma  s  5 1- Ma  5  5  0  )  /  Ma  5  5  0 

TD2  Di  5p2/ 62/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

rl  [61/50, 0,  LI/ 180,  0,  TDl/ 9,  0,  WHl/ 2,  2] 

r2  [  6  2  /  2  5,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

r2rl  [  6  2  /  2  5,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Heavel  =  seakeepi  nq(  '  FRMlh,  mat '  ,  rl); 

Rolll  =  5  e  a  ke  e  p  I  n  g  (  '  F  RMl  r ,  ma  t '  ,  rl); 

Fleave2  =  s  e  a  ke  e  p  i  n  g  (  '  F  RM2  h ,  ma  t '  ,  r2); 

Pitch2  =  5 e a ke e p i  n g (  '  F RM2 p ,  ma t '  ,  r2ri); 

$ 

=^PARAMETER 
$  n  a  me  value 

6rol  I  1  61>^Rol  11/57,3 

Lpitch2  L2>^Pi  tch2/ 57,  3 

znpr  ab5(6rolll)+ab5(Lpitch2)+ab5(Fleavel)+ab5(Fleave2) 

$ 

>^06]  VAL 

$  objective-function-value 
z  hpr 
$ 

>^END  06] 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  system  of  systems  level  optimization  - 
SOS  Constraint,  m 
6EGI  N  CONSTRAI  NT 

$ 

>^SYSTEM  OPTI  Ml  ZATI  ON 
$  s  y  s  t  e  m-  n  a  me 
HLSl 

$ 

>^SYSTEM  OPTI  Ml  ZATI  ON 
$  s  y  s  t  e  m-  n  a  me 
TS2 

$ 

>^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  SOS 

44441444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

$  Define  System  FILSl 
6EGI  N  SYSTEM 
$  s  y  s  t  e  m-  n  a  me 
HLSl 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  top  level  optimization 
6EGi  N  TOP 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  top  level  optimization  -  HL  S  1_  To  p_  Ma  i  n ,  m 
6EGI  N  MAI  N 

$ 

HI  NKI  NGVARI  A6LE 


$  n  a  me 

1  ower -  bounds 

upper  -  bounds 

61 

0,  9 

1,  1 

LI 

0,  9 

1,  1 

1 1 

0,  8 

1,  2 
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t  2 

0,  8 

1,  2 

1 3 

0,  8 

1,  2 

t  4 

0,  8 

1,  2 

t  5 

0,  8 

1,  2 

$ 

=^OPTI  Ml  ZATI  ON  OPTI  ON 
$  opt i  on  value 

Di  5  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

DIffMInChange  0,005 

OlffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

“^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  top  level  optimization  -  FILSl_Top_Obj  ect  I  ve,  m 
=^BEGI  N  OB] 

$ 

>^PARAMETER 
$  n  a  me  value 

Dh  12,0 

$ 

>^OBj  VAL 

$  objective-function-value 
(  Ll>^Bl=^Dh/  100,  0)  "0,  8/  2  6  7,  1 

$ 

=^END  OBj 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  top  level  optimization  -  HL  S  1_  To  p_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

>^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  displine-number 
Dl 

$ 

=^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  discipline-number 
D2 

$ 

>^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  TOP 

4444^444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

=^BEGI  N  Di  SCI  PLI  NE 
$  discTpline-number 
Dl 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  Dl  -  HL  S  1_  Dl_  Ma  I  n ,  m 
BEGI  N  MAI  N 

$ 


Hi  NKi  NGVARI  ABLE 
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$  name  lower-bounds  upper-bounds 
B1  0,9  1,1 

LI  0,9  1,1 

1 1  1,0  1,0 

1 2  1,0  1,0 

1 3  1,0  1,0 

1 4  1,0  1,0 

1 5  1,0  1,0 

$ 

=^OPTi  Ml  ZATI  ON  OPTI  ON 
$  opt  I  on  value 

Di  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

OlffMInChange  0,005 

OlffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  Dl  optimization  -  HL  S  1_  Dl_  Obj  e  c  1 1  v  e ,  m 
BEGI  N  OB] 

$ 

>^PARAMETER 
$  n  a  me  value 
WH 1  2,2 

Wei  ght  0  (  50,  0/  2M  0,  3  8  7  5  +0,  2  8  2  5  )  +1  2  M  0,  2  8  0  5  +0,  2  8  0  0  +0,  3  3  5  0  )  )  M80,  0 
Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 
DIspH  6  1  2  3  9,9 

Dlspl  DlspHMWelghtl-WelghtO)/WelghtO 

TDl  Dl  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

r  [Bl/50,0,  LI/ 180,  0,  TDl/ 9,  0,  WHl/ 2,  2] 

$ 

HTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Roll  =  s e a ke e p  I  n g  (  '  F RMl r ,  ma t '  ,  r) 

$ 

HB]  VAL 

$  obj  ective-functi  on-value 
Rol I / 0,  381 

$ 


HND  OBj 
$ 


$  Define  constraint  function  of  discipline  Dl  optimization  -  HL  S  1_  Dl_  Co  n  s  t  r  a  I  n  t ,  m 
HEGI  N  CONSTRAI  NT 
$ 

HARAMETER 
$  n  a  me  value 

WHl  2,2 

Wei  ght  0  (  5 0 ,  0/  2 M  0 ,  3  8  7  5  +0 ,  2  8  2  5  )  +1  2  M  0 ,  2  8  0  5  +0 ,  2  8  0  0  +0 ,  3  3  5  0  )  )  H 8 0 ,  0 
Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 
DIspH  6  1  2  3  9,9 

Dlspl  DlspHMWelghtl-WelghtO)/WelghtO 

TDl  Dl  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

r  [Bl/50,0,  LI/ 180,  0,  TDl/ 9,  0,  WHl/ 2,  2] 
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“^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Fleave  =  seakeepi  ng(  '  FRMlh,  mat '  ,  r) 

$ 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Pitch  =  5  e  a  ke  e  p  I  n  g  (  '  F  RMl  p ,  ma  t '  ,  r) 

$ 

=^CVAL 

$  I  nequal  I  ty-constrai  nt-functi  on-val  ue 
Fie  a  V  e  -  0 ,  17  0 
Pi  t  ch- 0,  155 

$ 

=^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  Dl  SCI  PLI  NE 


=^BEGI  N  Dl  SCI  PLI  NE 


$  discrpline-number 
D2 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  D2  -  FIL  S  1_  D2_  Ma  i  n ,  m 
BEGI  N  MAI  N 

$ 

$ 

HI  NKI  NGVARI  ABLE 

$  name  lower-bounds  upper-bounds 
B1  1,0  1,0 

LI  1,0  1,0 

t  1  0,8  1,2 

1 2  0,8  1,2 

1 3  0,8  1,2 

1 4  0,8  1,2 

1 5  0,8  1,2 

$ 

$ 

HPTI  Ml  ZATI  ON  OPTI  ON 
$  opt  I  on  value 

Dl  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

DiffMinChange  0,005 

DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 


$  Define  objective  function  of  discipline  D2  optimization  -  FIL  S  1_  D2_  Obj  e  c  t  i  v  e ,  m 
HEGI  N  OB] 

$ 

HARAMETER 
$  n  a  me  value 

wei  ghtl  (  Bl/ 2M  t  1-Ft  2)  -fUM  1 3-Ft  4-Ft  5)  )  HI 
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=^OBJ  VAL 

$  objective-function-value 
we  i  g  h  1  1  /  4  9  5  0 

$ 

=^END  OB] 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  D2  optimization  -  HL  S  1_  D2_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

>^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Secmod=SM(  t  1 ,  1 2 ,  1 3 ,  t  4 ,  1 5 ,  B 1 ) 

$ 

>^CVAL 

$  i  nequal  i  ty-constrai  nt-functi  on-val  ue 
2  2  7,  1  5  -  S  e  c  mo  d 


$  Define  System  TS2 
=^BEGI  N  SYSTEM 
$  s  y  s  t  e  m-  n  a  me 
TS2 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  top  level  optimization 
=^BEGI  N  TOP 
$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  top  level  optimization  -  TS  2_  To  p_  Ma  i  n ,  m 
BEGI  N  MAI  N 

$ 

HI  NKi  NGVARI  ABLE 


$  n  a  me 

1  0  we  r 

-  bounds 

upper  -  bounds 

B2 

0, 

90 

1,  10 

L2 

0, 

90 

1,  10 

El 

0, 

8 

1,  2 

E2 

0, 

8 

1,  2 

E3 

0, 

8 

1,  2 

r  hoi 

0, 

8 

1,  2 

r  ho2 

0, 

8 

1,  2 

rho3 

0, 

8 

1,  2 

HPTI  Ml  ZATI  ON 

OPTI  ON 

$  option 

value 

Di  s  pi  ay 

of  f 

Di  agnost i  cs 

of  f 

LargeScal  e 

of  f 

Ma  X  F  u  n  Ev  a  1  s 

50 

DiffMinChange  0,005 
DiffMaxChange  0,10 
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TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

>^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  top  level  optimization  -  TS2_Top_Obj  ect  i  ve,  m 
BEGI  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 

Dt  16,0 

$ 

>^OBj  VAL 

$  objective-function-value 
(  L2>^B2=^Dt/  100,  0)  "0,  8/  191,  25 

$ 

>^END  OBj 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  top  level  optimization  -  TS  2_  To  p_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

=^Di  SCI  PLi  NE  OPTi  Ml  ZATI  ON 
$  displine-number 
D1 

$ 

=^DI  SCI  PLi  NE  OPTI  Ml  ZATI  ON 


$  discipline-number 


44444444444444444444444444444444444444444444444444444444444444444444444444444444 

=^BEGI  N  Dl  SCI  PLI  NE 
$  discFpline-number 
Dl 


$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  Dl  -  TS2_  Dl_  Ma  i  n ,  m 
BEGI  N  MAI  N 

$ 

Hi  NKI  NGVARi  ABLE 

$  name  lower-bounds  upper-bounds 


B2  0,9  1,1 

L2  0,9  1,1 

r  hoi  1,0  1,0 

r  ho2  1,0  1,0 

r  ho3  1,0  1,0 


$ 

HPTi  Ml  ZATI  ON_OPTI  ON 


$  opt i  on  value 

Dl  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 


Page  7 


ca5e_5t  udy[  1] ,  t  xt 


OiffMinChange  0,005 
DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  D1  optimization  -  TS  2_  Dl_  Obj  e  c  t  i  v  e ,  m 
BEGI  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 
WH2  2,2 

A  1  3  6  9  8,  8  2=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  1  7=^B2H  2  /  2  5,  4  /  1  7  5,  0 

A3  3  8  4  8,  2  0=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  0=^  3  6  9  8,  8  2  +7  8  0  0,  0>^  5  4  4  9,  1  7  +7  8  0  0,  0>^  3  8  4  8,  2  0 

Massl  r  h  0 1  A1 +r  h  0  2  A2 +r  h  0  3  A3 

Di  spT  2  2  9  2  9,  1 

Di5p2  Di  5  pTM  Ma  5  5 1- Ma  5  5  0  )  /  Ma  5  5  0 

TD2  Di  5p2/ B2/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

r  [B2/25,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

$ 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Tpitch  =  5 e a ke e p i  n g (  '  F RM2 p ,  ma t '  ,  r) 

$ 

=^OBj  VAL 

$  objective-function-value 
Tpi  t  ch/ 0,  111 

$ 

=^END  OBj 
$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  DI  optimization  -  TS  2_  Dl_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

=^PARAMETER 
$  n  a  me  value 
WH2  2,2 

A  1  3  6  9  8,  8  2=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  1  1  *^2*12125.  4/  175,  0 

A3  3  8  4  8,  2  0>^  B2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  0=^  3  6  9  8,  8  2  +7  8  0  0,  0>^  5  4  4  9,  1  7  +7  8  0  0,  0>^  3  8  4  8,  2  0 

Massl  r  h  0 1  A1 +r  h  0  2  A2 +r  h  0  3  A3 

Di  spT  2  2  9  2  9,  1 

Disp2  Di  s  pTM  Ma  s  s  1- Ma  s  s  0)  /  Ma  s  s  0 

TD2  Di  sp2/ B2/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

r  [B  2  /  2  5,4  0,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

$ 

>^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Theave  =  seakeepi  ng( '  FRM2h,  mat '  ,  r) 

$ 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Troll  =  s e a ke e p  i  n  g  (  '  F RM2 r ,  ma t ' ,  r) 

$ 
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>^CVAL 

$  i  nequal  i  ty-constrai  nt-functi  on-val  ue 
Theave- 0,  343 
Trol I  -  0,  300 

$ 

=^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

>^END  Dl  SCI  PLi  NE 


44444444444444444444444444444444444444444444444444444444444444444444444444444444 

BEGI  N  Dl  SCI  PLI  NE 

$  discTpline-number 
D2 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  D2  -  TS2_  D2_  Ma  i  n ,  m 
BEGI  N  MAI  N 

$ 

HI  NKI  NGVARI  ABLE 

$  name  lower-bounds  upper-bounds 
B2  1,0  1,0 

L2  1,0  1,0 

El  0,8  1,2 

E2  0,8  1,2 

E3  0,8  1,2 

r  hoi  0,8  1,2 

r  ho2  0,8  1,2 

rho3  0,8  1,2 

$ 

HPTI  Ml  ZATI  ON  OPTI  ON 

$  opt i  on  value 

Dl  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 
DiffMinChange  0,005 

DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

HND  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  D2  optimization  -  TS  2_  D2_  Obj  e  c  t  i  v  e ,  m 
HEGI  N  OB] 

$ 

HARAMETER 


n  a  me 

value 

A1 

3  6  9  8,  8  2  H2H  2  /  2  5,  4  /  1  7  5,  0 

A2 

5  4  4  9,  1  7  H2H  2  /  2  5,  4  /  1  7  5,  0 

A3 

3  8  4  8,  2  0  H2H  2  /  2  5,  4  /  1  7  5,  0 

Wei  g h 1 2 

(  r  holHl-Fr  ho2H2-Fr  ho3H3)  H,  0  7  6  2 

$ 

HB]  VAL 

$  objective-function-value 
We  i  g  h  1  2  /  7 ,  7  2  4  4  e  4 

$ 

HND  OBj 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  D2  optimization  -  TS  2_  D2_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

=^PARAMETER 
$  n  a  me  value 

pres  0,  2=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

r  [El/3,  0el4,  r h o  1  /  7  8  0  0 ,  0 ,  E 2/ 3 ,  Oe  1 4 ,  r h o  2  /  7  8  0  0 ,  0 ,  E 3/ 3 ,  Oe  1 4 , 
rho3/ 7800,  0,  pres/ 0,  2] 

$ 

>^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

MaxMSl  =  MaxStress(  '  FRMl,  mat'  ,  r) ; 

MaxMS2  =  MaxStress(  '  FRM2,  mat'  ,  r) ; 

MaxMS3  =  MaxStress( '  FRM3,  mat'  ,  r) 

$ 

“^CVAL 

$  i  nequal  i  ty-constrai  nt-functi  on-val  ue 
MaxMSl- 250,  0e6 
MaxMS2- 250,  0e6 
MaxMS3- 250,  0e6 
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$  SS-MDO-U  analysis  with  Uncertainty 
$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

$  Define  system  of  systems  level  optimization 
BEGI  N  SOS 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  system  of  systems  level  optimization  -  SOS_Main,m 
BEGI  N  MAI  N 

$ 

=^DESI  GNVARI  ABLE 

$  name  initial-value  lower-bound  upper-bound  nondimensionalization-facor 
Standard-devi  ati  on 


Bl 

50,  0 

0 

8H0,  0 

1 

2H0,  0 

50,  0 

LI 

180,  0 

0 

8H80,  0 

1 

2H80,  0 

180,  0 

t  1 

0,  3875 

0 

7H,  3  8  7  5 

1 

3H,  3  8  7  5 

0,  3875 

t  2 

0,  2825 

0 

7H,  2  8  2  5 

1 

3H,  2  8  2  5 

0,  2825 

t  3 

0,  2805 

0 

7H,  2  8  0  5 

1 

3H,  2  8  0  5 

0,  2805 

t  4 

0,  2800 

0 

7H,  2  8  0  0 

1 

3H,  2  8  0  0 

0,  2800 

t  5 

0,  335 

0 

7H,  3  3  5  0 

1 

3H,  3  3  5  0 

0,  3350 

B2 

25,  4 

0 

8H5,  400 

1 

2H5,  400 

25, 400 

L2 

175,0 

0 

8H75,  00 

1 

2H75,  00 

175,00 

El 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

E2 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

E3 

3,  0el4 

0 

7H,  0el4 

1 

3H,  0el4 

3,  0el4 

r  hoi 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

r  ho2 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

rho3 

7800,  0 

0 

7H  8  0  0,  0 

1 

3H  8  0  0,  0 

7800,  0 

$ 

=^OPTi  Ml  ZATi  ON  OPT!  ON 
$  opt  i  on  value 

Di  s  pi  ay  iter 

Diagnostics  off 
LargeScale  off 
MaxFunEvals  300 
DiffMinChange  0,005 
DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  systemof  systems  level  optimization  - 
SOS  Obj  ect i  ve,  m 
BEGl  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 

WH 1  2,2 

Wei  ght  0  (  50,  0/  2M  0,  3  8  7  5  +0,  2  8  2  5  )  +1  2  M  0,  2  8  0  5  +0,  2  8  0  0  +0,  3  3  5  0  )  )  M80,  0 
Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 
DispH  6  1  2  3  9,9 

Displ  DispHMWeightl-WeightO)/WeightO 

TDl  Di  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

WH2  2,2 

A  1  3  6  9  8,  8  2  H2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  1  7  H2H  2  /  2  5,  4  /  1  7  5,  0 

A3  3  8  4  8,  2  0  H2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  OH  6  9  8,  8  2  +7  8  0  0,  0  H  4  4  9,  1  7  +7  8  0  0,  OH  8  4  8,  2  0 

Massl  r  h  0 1  OHl +r  h  0  2  0  H2 +r  h  0  3  0  H3 
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Di  spT  2  2  9  2  9,  1 

Di5p2  Di  5  pTM  Ma  s  5 1- Ma  5  5  0  )  /  Ma  5  5  0 

TD2  Di  5p2/ 62/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

rl  [61/50, 0,  LI/ 180,  0,  TDl/ 9,  0,  WHl/ 2,  2] 

r2  [  6  2  /  2  5,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

r2rl  [  6  2  /  2  5,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Heavel  =  seakeepi  nq(  '  FRMlh,  mat '  ,  rl); 

Rolll  =  5  e  a  ke  e  p  I  n  g  (  '  F  RMl  r ,  ma  t '  ,  rl); 

Fleave2  =  s  e  a  ke  e  p  i  n  g  (  '  F  RM2  h ,  ma  t '  ,  r2); 

Pitch2  =  5 e a ke e p i  n g (  '  F RM2 p ,  ma t '  ,  r2ri); 

$ 

=^PARAMETER 
$  n  a  me  value 

6rol  I  1  61>^Rol  11/57,3 

Lpitch2  L2>^Pi  tch2/ 57,  3 

znpr  ab5(6rolll)+ab5(Lpitch2)+ab5(Fleavel)+ab5(Fleave2) 

$ 

>^06]  VAL 

$  objective-function-value 
z  hpr 
$ 

>^END  06] 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  system  of  systems  level  optimization  - 
SOS  Constraint,  m 
6EGI  N  CONSTRAI  NT 

$ 

>^SYSTEM  OPTI  Ml  ZATI  ON 
$  s  y  s  t  e  m-  n  a  me 
HLSl 

$ 

>^SYSTEM  OPTI  Ml  ZATI  ON 
$  s  y  s  t  e  m-  n  a  me 
TS2 

$ 

>^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  SOS 

44441444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

$  Define  System  FILSl 
6EGI  N  SYSTEM 
$  s  y  s  t  e  m-  n  a  me 
HLSl 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  top  level  optimization 
6EGi  N  TOP 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  top  level  optimization  -  HL  S  1_  To  p_  Ma  i  n ,  m 
6EGI  N  MAI  N 

$ 

HI  NKI  NGVARI  A6LE 


$  n  a  me 

1  ower -  bounds 

upper  -  bounds 

61 

0,  9 

1,  1 

LI 

0,  9 

1,  1 

1 1 

0,  8 

1,  2 
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t  2 

0,  8 

1,  2 

1 3 

0,  8 

1,  2 

t  4 

0,  8 

1,  2 

t  5 

0,  8 

1,  2 

$ 

=^OPTI  Ml  ZATI  ON  OPTI  ON 
$  opt i  on  value 

Di  5  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

DIffMInChange  0,005 

OlffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

“^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  top  level  optimization  -  FILSl_Top_Obj  ect  I  ve,  m 
=^BEGI  N  OB] 

$ 

>^PARAMETER 
$  n  a  me  value 

Dh  12,0 

$ 

>^OBj  VAL 

$  objective-function-value 
(  Ll>^Bl=^Dh/  100,  0)  "0,  8/  2  6  7,  1 

$ 

=^END  OBj 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  top  level  optimization  -  HL  S  1_  To  p_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

>^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  displine-number 
Dl 

$ 

=^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  discipline-number 
D2 

$ 

>^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  TOP 

4444^444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

=^BEGI  N  Di  SCI  PLI  NE 
$  discTpline-number 
Dl 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  Dl  -  HL  S  1_  Dl_  Ma  I  n ,  m 
BEGI  N  MAI  N 

$ 


Hi  NKi  NGVARI  ABLE 
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$  name  lower-bounds  upper-bounds 
B1  0,9  1,1 

LI  0,9  1,1 

1 1  1,0  1,0 

1 2  1,0  1,0 

1 3  1,0  1,0 

1 4  1,0  1,0 

1 5  1,0  1,0 

$ 

=^OPTi  Ml  ZATI  ON  OPTI  ON 
$  opt  I  on  value 

Di  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

OlffMInChange  0,005 
OlffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^PARAMETER  U 

$  name  me  an  standard-deviation 

W1  2,2  0,2136 

$ 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  Dl  optimization  -  HL  S  1_  Dl_  Obj  e  c  1 1  v  e ,  m 
BEGI  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 
WH 1  2,2 

Wei  ght  0  (  50,  0/  2M  0,  3  8  7  5  +0,  2  8  2  5  )  +1  2  M  0,  2  8  0  5  +0,  2  8  0  0  +0,  3  3  5  0  )  )  M80,  0 
Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 
DIspH  6  1  2  3  9,9 

Dlspl  DlspHMWelghtl-WelghtO)/WelghtO 

TDl  Dl  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

r  [Bl/50,0,  LI/ 180,  0,  TDl/ 9,  0,  WHl/ 2,  2] 

$ 

HTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Roll  =  s e a ke e p  I  n g  (  '  F RMl r ,  ma t '  ,  r) 

$ 

HB]  VAL 

$  obj  ective-functi  on-value 
Rol I / 0,  381 

$ 

HND  OBj 
$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  Dl  optimization  -  HL  S  1_  Dl_  Co  n  s  t  r  a  I  n  t ,  m 
HEGI  N  CONSTRAI  NT 
$ 

HARAMETER 
$  n  a  me  value 
$  WHl  2,2 

Wei  ght  0  (  5 0 ,  0/  2 M  0 ,  3  8  7  5  +0 ,  2  8  2  5  )  +1  2  M  0 ,  2  8  0  5  +0 ,  2  8  0  0  +0 ,  3  3  5  0  )  )  H 8 0 ,  0 
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Wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5‘)  )  HI  ■ 

DispH  6  1  2  3  9,9 

Displ  DispHMWeightl-WeightOl/WeightO 

TDl  Di  spl/ Bl/ LI/  1  0  2  5,  0  /  0,  7  2  3  +9,  0 

$ 

HVAL  U 

$  i  n  e  q  u  a  I  i  t  y  -  c  0  n  5 1  r  a  i  n  t  -  f  u  n  c  t  i  0  n  -  V  a  I  u  e  t  a  r  g  e  t  -  r  e  I  i  a  b  i  I  i  t  y  -  i  n  d  e  x 
HLSlOlCons 1(  Bl,  LI,  TDl,  Wl)  2,06 


$ 

V  &p  u  1 

V  &p_  u  2 

V  5ip_  u  3 

v&p_u4 

v&p_u5  v&p_u6 

V  &p_  u  7 

V  &p_  u  8 

wr 

$ 

di  f  f  1 

di  f  f  2 

di  f  f  3 

di  f  f  4 

di  f f 5  di  f f 6 

di  f  f  7 

di  f  f  8 

0  ,  0  H  Wl 

* 

CVAL  U 

$ 

i  nequal  i  ty-constrai  nt-functi  on-val  ue 

target-rel i  abi 1 

i  t  y  -  i  n  d  e  X 

HLSlDlCon52(  Bl,  LI,  TDl,  Wl) 

2,  06 

$ 

V  &p  u  1 

V  &p_  u  2 

V  &p_  u  3 

v&p_u4 

v&p_u5  v5ip_u6 

V  &p_  u  7 

V  &p_  u  8 

wr 

$ 

di  f  f  1 

di  f  f  2 

di  f  f  3 

di  f  f  4 

di  f f 5  di  f f 6 

di  f  f  7 

di  f  f  8 

0 ,  0  H  W1 


HEGI  N  DI  SCI  PLI  NE 
$  discFpline-number 
D2 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  D2  -  HL  S  1_  D2_  Ma  i  n ,  m 
HEGI  N  MAI  N 
$ 

$ 

HI  NKI  NGVARI  ABLE 


$  n  a  me 

1  0  we  r  - 

bounds 

upper 

-  bounds 

Bl 

1,  0 

1, 

0 

LI 

1,  0 

1, 

0 

1 1 

0,  8 

1, 

2 

1 2 

0,  8 

1, 

2 

1 3 

0,  8 

1, 

2 

t  4 

0,  8 

1, 

2 

1 5 

$ 

0,  8 

1, 

2 

$ 

HPTI  Ml  ZATI  ON  OPTI  ON 

$  option 

value 

Di  5  pi  ay 

of  f 

Di  agnost i  cs 

of  f 

LargeScal  e 

of  f 

Ma  X  F  u  n  Ev  a  1  s 

50 

Di  ffMi  nChange 

0,  005 

Di  ffMaxChange 

0,  10 

To  1  Con 

1,  Oe-  2 

Tol  X 

1,  Oe-  2 

ca5e_5t  udy_u[  1] ,  t  xt 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  D2  optimization  -  HL  S  1_  D2_  Obj  e  c  t  i  v  e ,  m 
BEGI  N  OB] 

$ 

>^PARAMETER 
$  n  a  me  value 

wei  ght  1  (  Bl  /  2M  t  1+t  2)  +12M  1 3+t  4+t  5)  )  HI 

$ 

HB]  VAL 

$  objective-function-value 
weighti/49,50 

$ 

HND  OBj 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  D2  optimization  -  HL  S  1_  D2_  Co  n  s  t  r  a  i  n  t ,  m 
HEGI  N  CONSTRAI  NT 
$ 

HTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

Secmod=SM(  t  1 ,  1 2 ,  1 3 ,  t  4 ,  1 5 ,  B 1 ) 

$ 

HVAL 

$  i  nequal  i  ty-constrai  nt-functi  on-val  ue 
2  2  7,  1  5  -  S  e  c  mo  d 


$  Define  System  TS2 
HEGI  N  SYSTEM 
$  s  y  s  t  e  m-  n  a  me 
TS2 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  top  level  optimization 
HEGI  N  TOP 
$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  top  level  optimization  -  TS  2_  To  p_  Ma  i  n ,  m 
HEGI  N  MAI  N 
$ 

HI  NKI  NGVARI  ABLE 


$  n  a  me 

1  ower -  bounds 

upper  -  bounds 

B2 

0,  90 

1,  10 

L2 

0,  90 

1,  10 

El 

0,  8 

1,  2 

E2 

0,  8 

1,  2 

E3 

0,  8 

1,  2 

r  hoi 

0,  8 

1,  2 

r  ho2 

0,  8 

1,  2 

rho3 

0,  8 

1,  2 

$ 

HPTi  Ml  ZATI  ON_OPTI  ON 
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$  opt i  on  value 

Di  5  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

DIffMInChange  0,005 

DIffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^END  MAI  N 


$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  top  level  optimization  -  TS2_Top_Obj  ect  I  ve,  m 
BEGI  N  OB] 

$ 

=^PARAMETER 
$  n  a  me  value 

Dt  16,0 

$ 

=^OBj  VAL 

$  objective-functi  on-value 
(  L2=^B2=^Dt/  100,  0)  "0,  8/  191,  25 

$ 

=^END  OB] 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  top  level  optimization  -  TS  2_  To  p_  Co  n  s  t  r  a  I  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

=^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  displlne-number 
Dl 

$ 

=^DI  SCI  PLI  NE  OPTI  Ml  ZATI  ON 
$  discipline-number 
D2 

$ 

=^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

>^END  TOP 

4444^444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

=^BEGI  N  Dl  SCI  PLI  NE 
$  dlscFpllne-number 
Dl 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  Dl  -  TS2_  Dl_  Ma  I  n ,  m 
=^BEGI  N  MAI  N 
$ 

HI  NKI  NGVARI  ABLE 

$  name  lower-bounds  upper-bounds 


B2  0,9  1,1 

L2  0,9  1,1 

r  hoi  1,0  1,0 

r  ho2  1,0  1,0 

r  h  0  3  1,0  1,0 
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=^OPTI  Ml  ZATI  ON  OPTI  ON 
$  opt i  on  value 

Di  5  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 

DIffMInChange  0,005 
DIffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

=^PARAMETER  U 

$  name  me  an  standard-deviation 

W2  2,2  0,2136 

$ 

=^END  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  Dl  optimization  -  TS  2_  Dl_  Obj  e  c  1 1  v  e ,  m 
=^BEGI  N  OB] 

$ 

>^PARAMETER 
$  n  a  me  value 

WH2  2,2 

A  1  3  6  9  8,  8  2=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  1  1  *^2*12125.  4/  175,  0 

A3  3  8  4  8,  2  0=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  0>^  3  6  9  8,  8  2  ■F7  8  0  0,  0=^  5  4  4  9,  1  7  ■F7  8  0  0,  0=^  3  8  4  8,  2  0 

Massl  r  h  0 1 0”^  Al-Fr  h  0  2  0“^  A2 -Fr  h  0  3  0”^  A3 

Di  spT  2  2  9  2  9,  1 

Disp2  Di  s  pTM  Ma  s  s  1- Ma  s  s  0)  /  Ma  s  s  0 

TD2  Di  sp2/ B2/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  -F9,  0 

r  [B2/25,4,  L  2/ 1 7  5 ,  0 ,  TD2/ 9  ,  0 ,  WH2/ 2 ,  2  ] 

$ 

=^UTi  LI  TYPROGRAM  FUNCTi  ON 
$  function-call 

Tpitch  =  s e a ke e p i  n g (  '  F RM2 p ,  ma t '  ,  r) 

$ 

=^OBj  VAL 

$  objective-function-value 
Tpi  t  ch/ 0,  111 

$ 

=^END_OBj 


$ 


$  Define  constraint  function  of  discipline  Dl  optimization  -  TS  2_  Dl_  Co  n  s  t  r  a  i  n  t ,  m 
BEGI  N  CONSTRAI  NT 

$ 

=^PARAMETER 
$  n  a  me  value 
$  WH2  2,2 

A  1  3  6  9  8,  8  2>^  B2H  2  /  2  5,  4  /  1  7  5,  0 

A2  5  4  4  9,  Il*i2*l2l  2  5.  4  /  1  7  5,  0 

A3  3  8  4  8,  2  0=^  B2H  2  /  2  5,  4  /  1  7  5,  0 

MassO  7  8  0  0,  0=^  3  6  9  8,  8  2  ■F7  8  0  0,  0>^  5  4  4  9,  1  7  ■F7  8  0  0,  0>^  3  8  4  8,  2  0 

Massl  r  h  0 1 0”^  Al-Fr  h  0  2  0“^  A2 -Fr  h  0  3  0”^  A3 

Di  spT  2  2  9  2  9,  1 

Disp2  Di  s  pTM  Ma  s  s  1- Ma  s  s  0)  /  Ma  s  s  0 
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TD2  Di  5p2/ B2/ L  2  /  1  0  2  5,  0  /  0,  7  2  3  +9‘,  0 

$ 

=^CVAL  U 

$  i  n  e  q  u  a  I  i  t  y  -  c  0  n  5 1  r  a  i  n  t  -  f  u  n  c  t  i  0  n  -  V  a  I  u  e  t  a  r  g  e  t  -  r  e  I  i  a  b  i  I  i  t  y  -  i  n  d  e  x 
TS2DlCon5 1(  B2,  L2,  TD2,  W2)  2,06 

$  v5ip_ul  v&p_u2  v5ip_u3  v&p_u4  v&p_u5  v&p_u6  v&p_u7  v&p_u8 

W2  ■ 

$  di  ff 1  di  ff 2  di  ff 3  di  ff4  di  f f 5  di  ff 6  di  ff 7  di  ff 8 

0 ,  0 1  W2 

$ 

=^CVAL  U 

$  i  n  e  q  u  a  I  i  t  y  -  c  0  n  5 1  r  a  i  n  t  -  f  u  n  c  t  i  0  n  -  V  a  I  u  e  t  a  r  g  e  t  -  r  e  I  i  a  b  i  I  i  t  y  -  i  n  d  e  x 
TS2DlCon52(  B2,  L2,  TD2,  W2)  2,06 

$  v&p_ul  v&p_u2  v5ip_u3  v&p_u4  v&p_u5  v&p_u6  v&p_u7  v&p_u8 

W2  ■ 

$  di  ff 1  di  ff 2  di  ff 3  di  ff4  di  f f 5  di  ff 6  di  ff 7  di  ff 8 

0 ,  0 1  W2 

$ 

$ 

=^END  CONSTRAI  NT 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

=^END  DI  SCI  PLI  NE 

44441444444444444444444444444444444444444444444444444444444444444444444444444444 

$ 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 

=^BEGI  N  DI  SCI  PLI  NE 
$  discTpline-number 
D2 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  main  part  of  discipline  D2  -  TS2_  D2_  Ma  i  n ,  m 
BEGI  N  MAI  N 

$ 

HI  NKI  NGVARI  ABLE 

$  name  lower-bounds  upper-bounds 


B2  1,0  1,0 

L2  1,0  1,0 

El  0,8  1,2 

E2  0,8  1,2 

E3  0,8  1,2 

r  hoi  0,8  1,2 

r  ho2  0,8  1,2 

rho3  0,8  1,2 

$ 


HPTi  Ml  ZATI  ON  OPTI  ON 
$  opt i  on  value 

Di  s  pi  ay  off 

Diagnostics  off 

LargeScale  off 

MaxFunEvals  50 
DiffMinChange  0,005 

DiffMaxChange  0,10 
TolCon  l,0e-2 

Tol  X  1,  Oe-  2 

TolFun  l,0e-2 

$ 

HND  MAI  N 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  objective  function  of  discipline  D2  optimization  -  TS  2_  D2_  Obj  e  c  t  i  v  e ,  m 
HEGIN_OBj 
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=^PARAMETER 


n  a  me 

value 

A1 

3  6  9  8,  8  2>^B2H  2  /  2  5,  4  /  1  7  5,  0 

A2 

5  4  4  9,  1  7>^B2H  2  /  2  5,  4  /  1  7  5,  0 

A3 

3  8  4  8,  2  0'^B2H  2  /  2  5,  4  /  1  7  5,  0 

Wei  g h 1 2 

(  rhol'^Al-Prho2'^A2-Prho3'^A3)  =^  0,  0  7  6  2 

$ 

>^OBJ  VAL 

$  objective-function-value 
We  i  g  h  1  2  /  7 ,  7  2  4  4  e  4 

$ 

>^END  OB] 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$ 

$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

$  Define  constraint  function  of  discipline  D2  optimization  -  TS  2_  D2_  Co  n  s  t  r  a  i  n  t ,  m 
=^BEGI  N  CONSTRAI  NT 
$ 

>^PARAMETER 
$  n  a  me  value 

pres  0,  2>^  B2H  2  /  2  5,  4  /  1  7  5,  0 

r  [El/3,  0el4,  r h o  1  /  7  8  0  0 ,  0 ,  E 2/ 3 ,  Oe  1 4 ,  r h o  2  /  7  8  0  0 ,  0 ,  E 3/ 3 ,  Oe  1 4 , 
rho3/ 7800,  0,  pres/ 0,  2] 

$ 

=^UTI  LI  TYPROGRAM  FUNCTI  ON 
$  function-call 

MaxMSl  =  MaxStress(  '  FRMl,  mat'  ,  r) ; 

MaxMS2  =  MaxStress(  '  FRM2,  mat'  ,  r) ; 

MaxMS3  =  MaxStress( '  FRM3,  mat'  ,  r) 

$ 

>^CVAL 

$  i  nequal  i  ty-constrai  nt-functi  on-val  ue 
MaxMSl- 250,  0e6 
MaxMS2- 250,  0e6 
MaxMS3- 250,  0e6 
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